PIPS Lllffi EPPECIS Oil THE DYllilllC EESPOITSS OP IHS 

LIQUID 


A Thesis SuUirdtted 

In Partial Fulfilment of the Requirements 
For the Degree of 

MISTER OF TECMOIOG-Y 


vC'!?.. 

.,r 



pq P- PIP 

CHIKKAIAH PADMAITABHA 


to the 

Department of Mechanical Engineering 
lUDIM IlSTITUTE OF TE(SHNOLOGY, KAtTPUR 


|(\CLC- W 


V'"'- 

CM/ 1 f; 


May 1969 


Certified, that this work has keen carried 
out under my supervision and that the work has not 
"been submitted elsewhere for a degree. 



Dr. V. Srinivas 
Assistant Professor 


Dept, of Mechanical Engineering 



AGia;T0v3:iSDG3r.3nTS 


The author expresses his deep sense of gratitude 
to hr. V. Srinivas, "under whose able guidance this work 
has been completed. He would like to thank Dr. M.I.I.Oberai 
for his ever readiness to extend help, whenever it wa.s 
sought . 



1 


SYiropsis 


Pipe Line Effects On The Dynamic Response of the 

Liquid is a thesis submitted in partial fulfilment 
of the requirements for the degree of M. Tech. , 
by C. Padmanabha to the Department of Mechanical 
Engineering, Indian Institute of Technology, Kanpur 
in -the month of May, 1969 * 


This thesis analyses the effect of the longitudinal and 
transverse waves along a soft, elastic tube on the frequency 
response of pipe line flow. 


^\JChe exciting force for the longitudinal waves along thei^- 


tube is assmned to be given by »the shear drag oscillations on 




the walls of the tube'|.N>4he simplified tube momentum equation ' ' 

is solved for the above exciting force by considering the tube 

^ : 

as rigidly fixed at both end s'^ The exciting force is included 

in terms of Dirac delta function, the magnitude of which is deter-- 
mined by the exciting force at any point considered. Then the 
total effect of the exciting forces at'all points are included 


V,#' 

by integrating w.r.t.x from 0 to L'^^-dfter obtaining the longitu- 
dinal amplitude and the velocity of the tube, fluid momentum 
equation in the axial direction is solved for the velocity per- 
turbations by considering the velocity of the tube as a boundary 
condition at the inner radius of the tube. Perturbations in 
velocity and pressure response are added to the main velocity 
and pressure response, to obtain the total d 3 mamic response of 
the fluid flow. o/' 
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The exciting force for the transverse waves along the 
tube is assumed to be given by the force due to the pressure 
fluctuations on the walls of the tube. v::^he amplitude and the 
velocity of the tube in the transverse direct inn is obtained 
for the above exciting force as explained for the case of longi- 
tudinal waves along the tube'. The simplified momentum equation 
in the radial direction is solved for the radial perturbations 
velocity by considering the radial velocity of the tube as a 
boundary condition at the inner radius of the tube. Then the 
simplified continuity equation is solved to obtain the axial 
velocity perturbations by assuming zero velocity perturbation 
at the input end of the pipe line. With this as a boundary 
condition along the axis of the tube, the axial fluid momentum 
equation is solved to obtain axial velocity perturbations. Total 

response of the fluid is obtained by adding the main response to 

/ 

the perturbations response. 


^i^he 


results show^ that the effect of the longitudinal wave 
along the tube is very small as compared to effect of the trans- 
verse waves along the tube on the dynamic response of the fluid 
for a particular value of pipe line damping. The effect of the tu 

seemed to be increased Y/ith the increase of viscosity o-f the liqu 
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IRTRODUCTION 
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1.1 Physical Considerations 

Hydraulic transmission lines have been used in many 
industrial, research and military fields for high performance 
automatic control systems. The hydraulic servomechanism®, ope- 
rating at high pressure, are found more useful for high per- 
formance automatic control systems, especially in the air 
craft and rocket applications for the following reasons. 

1. Very high power amplications are obtained with 
such systems. 

2. Very high speed of response are possible. 

3. High power to weight ratio characteristics of 
these systems. 

The hydraulic liquids like water, oil etc. are used in 
measuring instruments like monometers, flow meters, transdu- 
cers^etc. Dynamic measurements made by these instruments are 
likely to have some error due to the dynamics of the liquid, 
elasticity of the interconnecting tubes etc. 

The effect of interconnecting lines may usually be 
neglected for' the low speed operations and measurements. But 
for high performance and high speed operations the effect of 
interconnecting lines may be predominant ’it is particularly 
so in the use of elastic tube applications’ such that the system 
performance may change significantly. The effect of interconnect 
ing lines may be important at all frequencies, when long or mo- 
derately long pipe lines are used. 
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If the fluid in these pipe lines were to he inTiscid 
and absolutely incompressible and the lines rigid, the sig- 
nal transmission through the fluid medium would ha-ye been 
instantaneous with the whole fluid medium in the pipe line 
acting as a single column. But, since no real fluid is in- 
compressible and inviscid, this is not possible in practice. 

Any disturbance of the fluid at any point gives rise to waves 
carrying these disturbances, and these waves travel in all 
directions. If the radius of the pipe line is small compared 
to the length of the tube, the effect of the disturbance tra- 
velling in the radial directions is very small since the radial 
waves are quickly damped aut because of high frequency radial 
fluid oscillations. These effects change the output response 
significantly as compared to the signal applied at the input 
end of the pipe line. 

The compressibility of the liquid together with its 
inertia gives rise to a finite velocity for transmission of 
small signals, called the phase velocity or the velocity of . 
wave, which causes certain time delay in signal transmission 
in addition to the dynamic behaviour of the output response. 

For small disturbances the phase velocity can be calculated 
from the linear wave equations for both undamped and damped 
sound waves. Knowing the phase velocity, the transient or 
the frequency response of liquid in the pipe line can be esti- 
mated by using the linear wave equations applicable to liquids 
excited by small disturbances. For large disturbances and 



3 


high pressure operations, the transient and freq[uency response 
are governed hy nonlinear simultaneous partial differential 
equations and, hence the analysis becomes more complicated. 

In this case phase velocity depends not only on the compressi- 
bility and the density, but also on the pressure level, the 
pressure disturbance amplitude, the rate at which pressure 
are varied, the flow velocity and viscous friction. 

Y/ith finite disturbance - signals applied to the liquid 
in a tube, waves may be transmitted in several modes and each 
mode has a unique wave velocity and damping characteristics. 

The major mode of wave transmission is due to the interaction 
of the inertia of the liquid and the capacitance effect of a 
compressible liquid in an elastic tube. Other wave modes that 
have been observed are the transverse, the longitudinal and the 
flexural' wave transmission in the elastic tubes. The trans- 
verse waves are noticeable in tubes having relatively low elas- 
tic modulus as compared to the compressible property of the 
fluid or in pipe lines subjected to extreme fluid pressure 
variations. The longitudinal waves are noticeable in tubes 
which are elastically soft and are induced either by the shear 
drag oscillations at the wall or by a sudden change in area 
of the pipe line at the ends. The flexural waves in pipe 
lines, which show oscillations in a direction perpendicular 
to the axis of the pipe line, depend on the type of line supports 
and are caused to some extent by the longitudinal oscillations 
at the end of a curved pipe line. 
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Ihe problem of the linear analysis can be analysed by 
several methods. A few most useful methods of analysis are 
the method of separation of variables, the Laplace transfor- 
mation method and steady state analysis. The Laplace transfor- 
mation can be used both for the transient analysis and the 
steady state analysis. The separation of variables can also 
be used for the transient case, but it is well suited for the 
steady state analysis. These methods analyse the response of 
the :fluid medium in pipe lines from the linearized equations. 

There are two methods for analysing the signal attenua- 
tion and dispersion effects caused by the shear stresses at 
the tube walls due to the viscosity of the fluid. The first 
method consists of solving the continuity and momentum equa- 
tions without neglecting. the radial velocity gradients in the 
shear stress terms contained in the momentum equation. The 
solution for the linear damped wave equation is obtained in 
terms of the modified Bessel function. The second method con- 
sists of replacing the shear stress terms in the momentum equa- 
tion by a term which is a product of the linearized resistance 
coefficient and the cross sectional average flow velocity. Bor 
a given tube and fluid, it is possible to find a relationship 
between the frequency and the linearized resistance coefficient 
using the linear theory for the oscillatory laminar fluid flow 
in pipe' line. In addition to the above effects, damping is 
Icnown to have a noticeable effect on the phase velocity and 
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hence on the resonant frequencies, linear damped wave equa- 
tions can he used to analyse the effect of the viscosity varia^ 
tions and of the linearized coefficient variations on the 
effective phase velocity. 

In soft elastic tubes, the gradients in axial flow velo- 
city perturbations which are caused by either radial velocities 
in the fluid due to the transverse oscillations of the tube 
walls or by longitudinal oscillations in the tube walls, intro- 
duce an additional effect on the dynamic flow of the fluid. 

These perturbations in the radial gradient of axial flow velocity 
introduce additional damping effects in the dynamic flow of fluids 
in tubes (Ref. 14). The radial gradients in axial flow velocity 
and hence the damping effects can be further influenced by the 
flexural mode of wave transmission in the pipe line. 

Ror the linear analysis of the problem of the dynamic 
response of the liquid in pipe lines, the properties of the 
liquid such as the compressibility and the viscosity are assumed 
to be constant. But these properties are known to vary consi- 
derably with temperature and pressure. In a system where the 
flow velocities and disturbances are resonably small, the effect 
of the temperature on the variation of viscosity and compressi- 
bility is not usually of much concern in the analysis of the dy- 
namic response of liquid in a pipe line because the liquid and 
the line has comparatively high specific heats which give rise 
to very small temperature changes in the system due to negligi- 
ble energy changes. 
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1.2 Earlier Work : 

Origin of the study of sound waves in long tubes has 

1 * 

been studied at the time of Kirchoff . later the transmission 
of sound in small tubes for various end conditions was studied 
by R^eigh . Since then, a great deal of work has been carried 
out on transmission of the sound waves in tubes. 

The purpose of early investigation on unsteady fluid 
flow in straight pipes was to analyse the surging phenomena, 
or water hammer, in power plants employing large diameter con- 
duits and to determine the velocity of sound in the fluid. The 
classical work in this direction was carried out by Joukowsky 
and llleivi'^ around 19OO1 

Wood-^ solved the water hammer problem using a damped wave 

equation consisting- of the linearized resistance coefficient, 

6 

assuming steady Poiseuille flow. Rich has given solution for 

damped (using linearized resistance coefficient) and undamped 

wave equations for several cases using the method of Laplace 

transformation. This method is more suitable for estimation of 

7 

both the transient and the frequency response. Iberall has 
given an interesting analysis for the behaviour of fluid in the 
transmission line with a pressure sensing instrument at the 
load end. Iberall then modified the analysis to add the effects 
of fluid compressibility, finite pressure excess, inertia of the 
fluid medium, finite length of tubing and nonlinear effects in 

* Superscript numbers refer to references. 
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Q 

the theoretical analysis. Nichols has analysed the problem 
for a semi- infinite pneumatic transmission lines for frequency 
response assuming the properties of fluid to be linear. He 
made use of the term called shin effect for the analysis of 
the frequency dependent viscous damping effects, to estimate 
the attenuation and dispersion effects caused by the viscosity 
of the fluid. This analysis estimates the frequency dependent 
viscous damping effects by neglecting the compressibility effects 
of the transmission medium in the pneumatic lines. 

Thomson^ has given an analysis for the existence of the 

several modes of wave transmission in liquids in a semi- infinite 

tube by considering the two dimensional undamped wave equation. 

His experiment gives no confirmation of the existence of these 

waves and he shows that below the cut off frequency of the first 

and the second mode these waves do not exist whereas slightly 

above the cut off frequencies the phase velocity for these 

10 

waves becomes infinite. Morgan and Kiely have given a method 
to analyse the effect of viscosity of liquids and the internal 

damping in a thin walled semi- infinite elastic tube on the 

11 

phase velocity. Gerlach and Parker have given an analysis 
to investigate the symmetric modes of wave propagation for a 
viscous, compressible liquid in a cylindrical conduit . Boun- 
dary conditions for both rigid and elastic walls are imposed 
and the resulting characteristic equation is solved for the 
spatial attenuation factor and phase velocity for several modes. 
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12 

Brovm has giTen an analysis for the dynamic response 
of pipe lines in terms of Laplace operator s and the modified 
Bessel’s functions. As this solution is obtained in terms of 
Laplace operator, it is valid for both frequency and transient 
response analysis. The computed results of the above analysis 
used approximate relationships for the modified Bessel’s func- 
tion, hence the results presented for the phase velocity against 
the . nondimensional frequency are not accurate for : noridimen- 
sional frequencies smaller than ,9. 

1'5 

D’ Souza and Oldenburger have given a method of esti- 
mating the frequency response of a long pipe line including the 
end effects. They have considered the longitudinal mode of pipe 
line oscillations as caused by the pipe line elasticity and the 
physical characteristics of the load block such as its inertia 
and damping. These longitudinal oscillations at the load block 
caused by the area difference between the area of the pipe line 
and the area of the load block on which fluid pressure is acting. 
Hence this problem is set up in a manner that holds good only 
if there is an area difference betv/een the area of the pipe line 
and the effective area of a freely supported load system, and 
if the pipe line is straight. 

14 . 

Srinivas ^ has developed the relationship betv^een the 
linearized resistance coefficient and frequency response of 
flow velocity.' It has been shown that the t?ro methods using 
the linearized resistance coefficient and the linearized vis- 
cous flow analysis for flow are the same when the frequency 
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variations in the linearized resistance coefficient as obtained 
from the developed relationship betv/een the tvro are considered. 
Attempt has "been made to include the other modes of wave transmi. 
ssion as caused by the longitudinal oscillations in the pipe 
line excited by the shear drag oscillations on the pipe wall 
and the transverse oscillations in the pipe line excited by 
pressure fluctuations in the fluid. An approximate method for 
the estimation of the additional resistance coefficient caused 
by the transverse vibrations in the pipe line has also been 
developed. 
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1*3 Object of this Investigation : 

Review of the earlier v/ork shows that much has been 
done on the theoretical and the experimental study of the 
dynamic response of the fluid medium in pipe line for a few 
restricted cases. But the effect of the wave propagation 
along the tube material on the dynamic response of the fluid 
has not been solved in detail. The effects of the tube for 

1 '^ 

the linear case have been analysed by D’ Souza and Oldenburger , 
assuming the exciting force is due to the area difference at 
the load end. Attempt has been made to include the effect of 
the wave propagation along the tube caused by shear drag osci- 
llations and internal pressure fluctuations in the tube by 
Srinivas”*^. 

This thesis makes use of the frequency response analysis 
for the dynamic response of the fluid flow. The effects of the 
longitudinal and transverse waves are included separately by 
assuming that the tube is rigidly fixed at both ends. The 
exciting force is included in terms of Dirac delta function, 
the magnitude of which is determined by the exciting force at 
any point considered. Then the fluid momentum and continuity 
equations are solved for the velocity perturbations and pressure 
perturbations. To get the total response, the velocity pertur- 
bations and pressure perturbations are added to the fluid velo- 
city and pressure response. The constants are evaluated by 
assuming suitable end conditions. 
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Chapter II giires the simplified fluid momentum and 
continuity equations.. Appendix A gives the detailed study of 
the order of magnitude of terms in the fluid momentum and 
continuity equations. Chapter III gives the frequency res- 
ponse analysis for the d3m.amic flow of the fluid, as given 
in the previous litera,tures Chapter IV gives the ana- 

lysis to include the effect of the longitudinal waves along 
the tube on the dynamic response of the fluid. Chapter V 
gives the analysis to include the effect of the transverse 
waves along the tube on the dynamic response of the fluid. 
Chapter VI gives the detailed discussion of the computed 
results. 



CHiPTER II 


ASSTOIPTIOES MD BASIC EQUATIONS EOR THE EIUID ELOW 



2.1 Differential Equations for the Fluid Flow: 


As this analysis is dealing with the pipe flow, it 
is convenient to use the differential equations in cylindri- 
cal coordinates, let x he the axis along the axis of the 
tube. By assiming axial symmetry of the flow the deriva- 
tives and the velocity component in the circumferential 
direction can he neglected. Assuming the origin lies on 
the axis of the tube 


continuity equation is given by 
Jl ^ ^ j. ^ j. ^ 4. Z ^ 

f 9t ¥r r 9x i* 9x 5" 9x 

1 5 

Equation of motion is given by 


0 


( 2 . 1 ) 


9v ^ 
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9v 
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^ 9x 

+ u 
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, „ 9u 

Ft 


+ u 


1 As 

f 9x 


+ V 


9^v . 9^v J_ ^ 

ax2 ar2 




9r" 


9^u 

9x^ 


J _ u^ 

r 9r ^2 


(2.3) 


The isentropic equation for the liquid is assumed to 
be of the form 



(2.4) 


If the pipe line is very long compared to the inside 
diameter and if the pipe line is quite rigid, the radial velo 
city and its - gradients and the pressure in the radial 
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direction are damped out quickly. Hence, the radial irelocity 
and the momentum equation in the radial direction can he 
neglected. 1 detailed discussion of the order of magnitude 
of terms is given in Appendix A. Simplified form of the 
continuity and momentum equations are written as 


J. 4- 

s 9t 6x 


= 0 


(2.5) 


— +- V L ■— p + T 6TJ" ° 


6 1 


dx 


f - - 6r' 

Combining Equations (2.4) and (2.5) yields the 
following equation 


( 2 . 6 ) 


1 |£ 
B H 


+ = 0 
B ax ^ Sx ^ 


(2.7) 


Replacing the shear stress term in Equation (2.6) 
by an equivalent linearized resistance coefficient R, 
Equation (2.6) reduces to the form 


9 V 

at 


1 




0 


( 2 . 8 ) 


The above equation can be used to analyse the problem 
for the linearized resistance coefficient approach. 

The isentropic wave propagation velocity in a compre- 
ssible fluid for small perturbations is given by the equation 


= f dp/df 

Erom the equation of state of liquid_, 

dp/d? = B/p = Cq^ (2.9) 
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2.2 Effect lire Bulk Modulus : 


Due to the law compressihility of the liquid, the 
elastic property of the pipe line will have some influence 
on the effectiTe bulk modulus which in turn, affects the 
overall lesponse of the system. The elasticity of the tube 
material causes a reduction in effective bulk modulus which 
reduces the wave propagation of the liquid. A method of 
estimating the effective bulk modulus was given by Srinivas 
The method estimates the effective bulk modulus of liquid 
in thick walled pipe line when the variations in the fluid 
tube dilation are small as compared to its radius. The 
details of the above method is given in Appendix C, and the 
effective bulk modulus is approximated by the equation 


2R,- 


[(l-<r-) + (l+>') 


2 

i 


B’ 


E (R„‘ 


Rib 


14 

* 


± 

B 


( 2 . 10 ) 
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2.3 Assumptions : 

The following assumptions are made for the dynamic 
flomr of viscous liquid in a pipe line in order to facilitate 
a linear analysis of the general wave equations and the 
equation of motion. 

1 . Fluid flow is axially symmetrical and the end 
effects are relatively small for a long pipe 
line . 

2. The flow is laminar with constant fluid viscosity 
and hulk modulus. 

5. The hody forces are negligible. 

. The fluid compressibility is relatively small 
compared to the fluid oscillations. 
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CHiPTER III 


MilYSIS POR THE LINEAR CASE 
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3<,1 Frequency Response Analysis : 

From Equations (2.5) and (2.6) the simplified continuity 
and momentum equation for the linear case is given "by 


1 _ 

f 



+ 


dv 

ex 


= 0 


(5<.1) 


9 V 1 
ct dx 


dx 



dv 

dr 


0 


(3.2) 


By combining Equations (2.5) and (3.1) the following 
equation is obtained. 


1 9y 

B St dx 


0 


(5.3) 


The separation of variable method uses the substitution 
joofor d/dt to separate the time derivative and make it suitable 
for steady state response, analysis. The above substitution 
reduces the partial differential equation into an ordinary diffe- 
rential equation called the Helmoltz equation. 

It is assumed that the solution of Equations (3.2) and 
(3.3) be of the form 


p(x, t) = 


(3.4a) 

v(x,r,t) = 

V^(r) e^ 

(3,4b) 

u(x,r,t) = 

U3^(r) e^^ 

(3.4c) 


Substituting Equation (3,4) in Equation (3.2), Equation 
(3.2) reduces to the form 
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3^2 r ^ 


Y 

1) X 


r 


dP 

X 

dx 


0 


Substituting T 


. dP 

Y + — 1 — — — 

X jo)? dx 


(3.5) 

(3.6) 


Equation (5.5) reduces to the form 


^ 1 
ar^ ^ 


^ _ iw T = 0 

dr ^ 


(3.7) 


The above equation has only one solution which is 
finite at r = 0, and this is given by 

T = h(x, DW) iQ(Kr) (5.8) 

where is the modified Bessel's function of order 
zero 3j\di 


K = nwr~ 

Combining Equations (5.6) and (5.8), the following 


equation is obtained. 

= h(x, Dcn ) (Kr) 


1 (3P 

1 X 


X 


o 


3 oif dx 


(3.9) 


Boundary conditions are 


at 


r = R. 


0 


\(r = Rj_, X, gcb ) 

av /8r = 0 

IK) 


0 


(3.10a) 

(3.10b) 


Using the first boundary condition in Equation (3«9) 


yields 


d P. 


X 


dx 


gcqf h(x, geo ) I^(K R^) 


(3.11) 
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Uow 1 combining Equations (3. It) and (3.9) yield' 


= h{^, 3u)) [l^ (Er) - I^(K R^) ] (3.12) 

Equation (3.12) shows that the amplitude of axial flo?/ 
velocity has a maximum value at the centre and decreases to 
zero at the walls of the tube. The average cross sectional 
axial flow velocity and pressure can be defined by the follow- 
ing equations. 




/ ^ rdr 
0 
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uRi^ 


/ ^ r dr 
0 


Erom Equations (3.12) and (3.13) 


(3.13) 


(3.U) 


h(x, 3 cn) 


(3.15) 


where K, 


K R, 


( K R^) - (KR^) 


As it is assumed that the pressure is a function of 
X only, the average cross sectional pressure can be written as 

^x = \ 

Substitution of Equation (3.4) in Equation (3.3) yields, 


LI2 

dx 


i]. y. p 

B -^x 


(3.17) 


The average axial flow velocity is used in this equation 
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because of the assumption of axial symmetry and negligible 
radial Telocity, Differentiation of Equation (3.17) yields 


^ JB ^ 


(3.18) 


Substituting Equations (3.11) and (3*15) in Equation 
(3.18) the following equation is obtained. 


Tdiere 


h(x, jco ) 
dx2 


h(x, jcD) = 0 


I,(K R.) 
* K. 


(3.19) 


tion 


Solving Equation (3.19) for h(x, ooO ) yields the solu- 


h(x, joe) = G-^ cosh(f x)+ sinh(f x) (3.20) 

Substituting Equation (3.20) in Equation (3.15) the 
average emss sectional axial velocity component is given by 

^ 0^ cosh(J' x)+ sinh (f x)J (3.21) 

Erom Equation (3.17) pressure amplitude is given by 


p = 

X CO 


dx 


( 3 . 22 ) 


Substitution of Equation (3.21) in Equation (3.22) the 
pressure component is given by 


13 : 

X ' , 


sinh(r xj+ cosh(^if x}J (3.23) 


By neglecting the effect of tube and using Equations 
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(3*21) and (3 •23), and using the follo?/lng boundary conditions 




Y 


( X 


P ( X 

X ^ 


\ 


0 , 

0, 30 )) 

I>, D<^) 
Ijj 3 o3) 


= V. 




the constants and are obtained as 


G 


1 


¥ 

5 


H 


3 P 


1 


(3.24) 


(5.25) 


"2 ' 9 0^ (K R. ) 

“'0 0 1 

Then from Equations (3.21), (3.23), (3.24) and (3.25) 
the following equations for the pressure and the velocity 
response are obtained. 


¥ 

_c 

P 


cosh(ir Ii) - — ^ 


D 


(5.26) 


¥ 


jr 

cosh {( Jj)- T— D sinh(f Lj 


(5.27) 


where P 


Pj’c 


Dividing Equation (3.27) by Equation (3.26) and re- 
arranging the teims yields the following equation. 


P^ cosh (r l) + ^L/D sinHf l) 
¥^ ^ Pj^ 'D sinh(f Li)+ ¥^ cosh(f l) 


(3.28) 
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By substituting the above equation in Equations (3.26) 
and (3.27) and by further simplification, the frequency 
response equations for the pressure and the velocity are 
obtained as 

P 1 

—ii = (3.29) 

^o coBh(r L) - 5 py c^ sinh('r l) 


Y. 


1 


CO 


sh(r i) 




(5.30) 


Where is the load impedence which depends upon 

the condition at the load end of the pipe line. 1 method of 
estimating the load impedence for a load system with mercury 
column oscillating in a vertical U-tube is given in Appendix E, 
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3<.2 The Effect of Tube on the Dynamic Response of the Fluid : 

In addition to the change in effectiTe bulk modulus 
of the liquid, the wave transmission along the material of 
the pipe line influence the dynamic response of the fluid 
media. The wave transmission along the pipe line usually 
takes place in several modes. In the following chapters only 
the axial symmetrical modes are considered for the analysis. 
They are the longitudinal and the transverse modes of the 
pipe line vibrations. 



GHiPTER IV 


IKE EEEECT OE THE LOHGITUDIHil WAVES IN THE 
TUBE ON THE DYITMIC RESPONSE OE THE LIQUID 
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'The Effect of the longitudinal waves in the tube on the 
Dynamic r e s ponse of the liqui d : 


given by 


The tube momenttun equation in the axial direction is 
17 




G- r If ^ ^ damping 

6r 

where z is the longitudinal oscillation amplitude 
y is the transverse oscillation amplitude 


(4.1) 


The above equation is obtained by considering the 

compressibility of the tube material. The detailed study of 

Equation (4.1) is given in Appendix 1 for the frequency 

response analysis. The resulting wave equation is linearized 

by neglecting the terms of smaller order of magnitude and the 

18 

damping effect is included by adding viscovelastic damping. 
The resulting linearized damped wave equation is given by 




a^z 

at^ 


= CA+ 20) + (K+ 4->i) ^ 


ax^ 


at8x^ 


(4.2) 


Eor the frequency response analysis the solution of the 
equation is assumed to be of the form 


z(x,t) 


z e 

X 


D«*^t 


(4.3) 


Substituting Equation (4.3) in Equation (4.2) yields 


the following e quat i on . 
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=0 (4.4) 

r ^ = >, f 'zs\ ; El = K + H 

General solution of Eq_uation (4.4) without the forcing 
f^mction is given by 

= G^ cosy X + sin y x (4.5) 

Assuming the tube is fixed rigidly at both ends, the 
end conditions are 

at X = 0 2 _ = 0 (4.6a) 

X = L z = 0 ■ (4.6b) 

JL 

Substituting (4.6) in Equation (4.5) yields 

G^ = 0 and siny L = 0 (4.7) 

Eor nontrivial solution ^ 0 

Hence sin y I = 0 

or Y = ^ where n = 1, 2, (4.8) 

Substituting Equations (4.7) and (4.8) in Equation 
(4.5), Equation (4.5) reduces to the form 

Now considering only the excitations caused by the 
shear drag oscillations, the distribution of the shear forces 
on the tube wall can be estimated as given below. 
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The axial flow velocity distribution for viscous flow 
is given by Equation (5.12) as 

= h(x, 3 W) [l^CKr) - (KR^) ] (4.10) 

The shear stress on the tube wall is given by the 
equation 

8 V 

=Kh(x, 2 ^) K I.(ER ) (4.11) 

r=R^ ^ ^ 


let dx be an elemental length of the tube at any point 
x^, Heglecting the forces acting at all the points except at 
X = x^, the forcing function is given by 

8 V I 

Forcing function = 6(x - x ) dx 27i:R. K- (4.12) 

° 1 [r=R^ 

Where 6(x-x^) is Dirac delta function having the 
property 


6 (x - x^) =1 ■ when x = x^ 

6 (x - x^) = 0 when x 7 ^ x^ 

Row Equation (4.2) is multiplied by the area of cross 
section of the tube material and elemental length dx. Then 
this equation is equated to the forcing function (4.12). After 
rearranging the resulting equation the following equation is 



6(x-z^)2ti:R^K- 



9v t 

It I r=R^ 


(4.15) 
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where = area of cross section of the tuhe material 

c 

lahing Laplace transform of Equation (4.13) yields 


p — 

s z 


- CA+ 2G) ^ - (K + In) . ^ 

dr. ^ dx^ 


^6(x - x^) 2ii:R. f 


0 V 


r=R,- 


(4.14) 


By substituting s = jo) for the steady state response, 
Equation (4.14) reduces to the foim 

av 

6(x-x^) 2TtR.K-^ 


Now : substituting Equations (4.9) and (4. It) in 
Equation (4.15) » and ' rearranging the resulting equation 
yield 

S (.-fj + ) sin = 

h(x, 3 co) 6(x - x^) (4.16) 

where = (>+ 2G) + 3 co(V\+ 

2% R. H-K I. (m.) 

ri — ■ ■ X ■ i X 


a^s. 


3CJ 


dx" 




Now by the theory of Eourier series and by the property 
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of the Dirac delta fimction 


/ f(x) 6(x - X ) dx = f(x ) 

0 ° ° 

G-^ is obtained from Equation (4.16) as 


= 


20^ h (x^, 3co) sin (nicx^/L) 
I + a^n^it^/L^) 


(4.17) 


Substituting Equation (4.17) in Equation (4.9) the 
longitudinal amplitude of the tube due to the forcing func- 


tion at X ■= Xq is given by 


nmx 


z 


2C| h(xQ, 3(u>) 


sin 


° sin S# 


XX. 




+ 


2 2 2 
on % 


(4.18) 


Longitudinal oscillations of the tube due to all the 
forcing function acting on the tube from 0 to L is given by 

Erom Equations (4.18) and (4.19), z is given by 


where G, 


/n( Xy Itk ^ a 

^2 ^ ^n L 

4-Tt^ I|(KEj_) 

L^ A_ 


( 4 . 20 ) 


n 


n [ G| (1 - (-1)^ cosh(f L))-H^(-1)^ sinh(f L) 

(f^ + + a^n^Tz^/Jj^) 

( 4 . 21 ) 

Ihese details are given in Appendix B« The axial velocity 
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perturbations can be estimated by considering unit length 
of the tube as given below. 


Velocity amplitude of the tube = dz /dt = (4.22) 

-K. X , 

Let v’ be the axial velocity perturbations of the 
fluid due to the longitudinal oscillations of the tube. Then 
the fluid momentum equation for axial velocity perturbations 
is given by Equation (3.5) as 


5^v’ . ov’ 

X _1 X 

.2 r 




3.r 


__ 


dp’ 


X 


p- dx 


= 0 


(4.23) 


where p’ is the pressure perturbations. 


By going through the steps indicated in Equations 
(3.6), (3.7) and (3.8) the solution of Equation (4.23) is 
given by 

ap’ 


vq = g(x, ju.) q(Kr) - 


(4.24) 


at 


Boundary conditions are 
6v* 

^ _ o X = 0 


9r 


r = R. v' = jcuz from Equation (4.22) (4.25) 

XX X 

Now using the boundary condition at r = Rj_ in Equation 
(4.24), the following equation is obtained. 


dp’ 


X 


dx 


3^f t gCx> Iq (KR^) - 3coz^ i (4.26) 


Combining Equations (4.26) and (4.24) yields the follow- 


ing equation 



2 y 


-^' 3 , = g(x, DCO) [ I^CKr) - Io(KR.) ] + j U^2l] 

Average velocity peirturbations is obtained by the 
equation 


V 


2 % 


X 


or 


0 


= K. 


R. 

T* r d r 

£L r\ X 


(4.28) 


x -2 g(x, 3^4 + 

Ror the velocity perturbations Equation (3.18) can be 
written as 


dp’ 


dx tn 


d^^’. 

dx^ 


(4.29) 


Substituting Equations ( 4 . 26 ), (4.28) and ( 4 . 20 ) in 
Equation (4.29) and rearranging the terms yield- the follow- 
ing equation. 


g(x, jo.) - ‘' 2 ^ 


2_2 


dx" 


K, 


CO n It 
- 2 ( 


. 2-2 


( 4 . 30 ) 


These details are given in Appendix B. 


Particular integral of the above i nonhomogeneous 
ordinary differential equation is given by 


, (H^ . c^/o 2) 

Oo 3co- t .2 ^ o ^ 


Sr. = - 




(f2 + nV/L^) 


■h- HTCX 

0 ^ sin j- 


( 4 . 31 ) 


Hence the general solution of the differential equa- 
tion (4.30) is given by 


g(x, jco) = cosh(fx)+ sinh(fx)+ g. 


P 


(4.32) 
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Substituting Equations (4.20) and (4.32) in Equation 
(4.28) the cross sectional average velocity perturbations 


is given by 


X “ "^4 


Kg cosher x)+ Kg sinh(fx)+ Gg 0 '^(r^+ ) 


■h o-iv, 

’"n — 

ir + 

L 


(4.33) 


Because of the reason that' a '-unit length of the;' tube is 
considered in Equation (4.22) the average axial velocity 
perturbations at any point x is given by 

X _ 

V = J V ’ dx 
0 ^ 


V' = G 


^ sinh(r x)+ (cosh(fx_)- 1) + 

(1-cos nxx/l) 

^ 2 

n(r2+2^) 

Cg ( f 2 +a>2/cQ^)L 3 


(4.34) 


where C. 


Erom Equation (3.17) pressure perturbations at any 
point X is given by 

T 3 . = a® ^ 

^ cjO dx 


Substitution of Equation (4.34) in the above equation 
yields the following equation. 



pit: 


P’ = Gc cosh r X + H(- sinh f x + 


CO I, 


sin nitx/L 
(f ^+n^-n; 


(4.35) 


To solve for the two constants and following 
assumptions are made 

at x=-L^= 0 or ^=0 

dx^ 

As there is no further effect of the tube, at the 
end of the tube, the 'assumption that there is no further 
change in pressure perturbations, is justified 

at x=li p*=P'j^ 

V' = ¥' 


In addition to the above assumptions it is assumed that 


^1 


Vt 


By substituting the above end conditions in Equations 
(4.34) and (4.35) the following equations are obtained 

.2 


G^ sinhfl + coshfl 




(-1)® n b 


n 


2 2 

(r^ + a-s- ) 


C4-.36) 





G-r 


Y. 


L ^ cosh(r L) - sinhcr l)]+H^ [ ^ ^ sinh(f Ij - 
cosli(ri)+ 


P£ 


n 




(4.37) 


The above two eq_uations can be written in the matrix 


form as 


sinh^f l) 


f cosher l) 


V. 


f 

T oosh-if l) - sinh(ir l) 


Y» -TSV 

. ^ sinh(f l)- cosh(f l)H 

L 


/ 


w 


(- 1 )“ n t 


n 


H. 


(f^+ n^TiVl*^) 


I 


2C 


n 


nCf^+n^u^/L^) 


By applying the Cramer’s rule and are obtained 
from the above matrix as 


% = [ 20l1 =Osll(fL)g.; 


n 


n(f^+n^7cVl'^) 


+ 


Y' 


1^ " 1 
(-1 n b 


( ^ sinh(fl)- cosh(fL) + 1) 


n 






(4.38) 
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^5 = 


• 2 C^^ f 2 sinli(rL) 2 _. 


n 


nCr^+n^-n:^/!^) 


^2 Y? 

— -p --' ( ^ cosh(Y‘L)- sinh(f L)) 


■^L 




2 " P. 

C-1)“nt^ . 


(f^ + n^%^/'L^) 


-1 / r pr V« . 

-73- pT - slnhjr 1) 


( 4 . 39 ) 


The frequency response of the velocity including the 
effect of the longitudinal waves along the pipe line can he 
obtained by adding Equations (3.22) and (4.34). 


Y = T + V’ 

X 


( 4 . 40 ) 


By substituting the constants and in Equation 
(4.34) and then substituting Equations ( 4 . 34 ) and ( 3 . 22 ) 
in Equation (4.40) yield 

r r V 

"V = Gi L ”■"2 — sinh(f x) ( sinh(fl)- cosh(rij)n )S 

Ti Tt 


+ 2 r ^ cosh(r 1) sinh(f x) Sp - ^ sinh(fL)* 

C 

(cosh(f x) - 1) Sp - - (cosh(r 3^1 )(^ cosKfl) 

^ ^1-cos n%x/j^ \ 

- sinh(f l))S^ ] / 2 — 




cosh(f 1)) + Kp cosh^jf x) 


( continued on the next page . . . . ) 
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wlie re 


n , 2 y . 

+ sinh(fx)(^ ^ sinh(ir 4 -cosii(ri>-l)S 5 


¥' 


-2 cosli(rLj sinh(f x) + 2C^^ sinh(r35(cosh(r3^1 )S 

C 

4 . __11_ (cosh(r:^1) cosli(f L]-sinfcL(f l))S^ J /Orj 


-C 11 2 


(1-coa ^ 


(-1)^ sinh(/’L)+ K„ sinhCTx) 


(4.41) 


n 


Further simplification of Equation (4.41) yields 


^ = «nx1 


°x1 Sl + Sj 


mrx 


"nxl 




(1-oos j 


9 


(4.42) 


[ 1-C-1) 




cos. 


n 


+ cosli(f x^ 


nTcx 


•«xlS3 - 


'^nx2 


- 011 " 


^1-COS 


(-1 )^ sinh(/' J^y 


'1 


a„ 


n 


sinja(J x ) 

'x1 


Q ^ L ( sin]i(f 1 ) sinh^r x) - cosii(f x) cosli(r E) 

Ii ' L . i 


cosli^f £)) - sin]i(r x) cosli(fl)+ sinli(f l) cosh^f x) + sinhir:^sin]:!] 


Cx 2 = 2C^^ f ^ (cosii(fi) sinhCfx}- sinh(f l) cosh(f x) + sinh(f I 
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(- 1 )^ coshCn)] 


S, = E 


n 


cosii(f l) 


*/ 3 - “ 


; = S — — 

5 S'. 


sinh^L) 


n 


n 


(- 1 )^ sinhrL) 


S. = Z 


^4 


>.?--■ 


n 


a = 
n 


C., = 


7 


■f 

^ siniilf l) 


pf%.* 


a;pt 


at X = 0 


Cqi + Cq2 S 2 

^nOI - 57 -^ ^ ^2 


a - -"01 S - ^2^4 

°n 02 - 


^ 01 "- ~ ® ’ ^02 ^ 


at 


X = It 

nl 1 


"l 1 "l 2 ^2 

"7 


+ 2 C11 $2 + K2 cosli(f LJ 


-C,, S. - S, 


^nL2 


^ 12 2 _ 2c^_j + £2 sinli(r l) 


°L 1 = ^ (ooshfL)- 1)] 


t 

1 


^12 " ^ ^ sinhlfl) 
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The frequency response of the pressure including the 
effect of the longitudinal wares along the pipe line can he 
obtained by adding Eq.uations (3.23) and (4.35). 


P = + p' 


(4.43; 


By substituting the constants G-^ and in Equation 
(4.35) and then substituting Equations (4.35) and (3.23) 
in Equation (4.43) yield 




|~ _L1 — cosh(r x) (^ ■ sinh(f 1} - cosh(f L)+1 ] 


Y' 


^P* UP 


+ 20^^ cosh(lf e) cosh(f x)S 2 - 2 f^sinhlfE)- 

0 "jif ^ "\r ^ 

sinh(fx) $2 - sinh(f x) (^ cosh(ri)-sinh(n^)S^ 


.-n n sin (nTix/l) ' -i 

E (l-(-1)'^ cosh(f])) + Cg sinh(/x) 

P -i 
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+H. 


0 .A 

cosh(lf x) (^ sinh^I^ - coshCfl^l )S 


•2Ci^ cosh(f L)cosh(r x} + 20^^ sinh(rE3 ^ 




sinh(f x) + - - - 2 ■^^3’ 3c) (^ cosh(fl) 


- sinh(fL)) S J / 


■n-Ts.. n sin nxx/E 

^ z — ^ ^ -(-D^snihfi) 


n 


+ Gg cosh(f x) 


(4.44) 


V kJ 
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Further simplification of Equation (4.44) yields 
^ = °nx3 ®1 ^ =nx4 


where 


n _ ^x3 °x4 ^2 ^ ,1Btc 

^nx3 " ^11 cnL 

^7 


n[ 1-(-1)^ coshfl] sin 


n-n:x 

1 


n 


-C , S, - C„, S, 


^nx4 


+ Gg coshlf x) 


x3 "3 "x4 ‘-'4 n nB Ti 2 


+ Cg sinh(f x^ 

’ n(-1 )^sinhCL) sin ^ 


'7 


11 coL 


a 


n 


*^x3 ~ (cosh({x) sinh(f L)- sinh(f:^cosl{fl)) 


+ coshlf x)- cosh(/'x)cosh(/lj + sinh(jrx3sinh^l)1 
2 ^■ 

x4 cU 

3^c I (ER ) 


2C..f^ jBf - 

— -^ — (cosh (rp) cosh(fx)- sinh(r Iij sinh(rx)) 


0 


8 


at X = 0 


"n05 


°03 ^ Co4 Sp 

O 7 


c., s, - e., s. 

Q _02_3 04_4 0 

n04 n. 8 

^^7 



38 


at 


”''2 ” ^ ^ I) + 1 - ijl 


2C^^ f2 pf 


CO 


COSllCf l) 


"03 ■ 

=04 = 

X = L 


Ct, S. + C., S. , 

\l5 = + Cg slnn(f l) 

^7 


C, 


nL4 


%3 


-Ql3 ^3 ~ ^4 ^4 


'7 
.2 _•■ 


+ Gg coshCfl) 

G.p" jBf 

— (coshCn)- 1) 


C. 


L4 


DBf 

__ 


Using Equations (4.45) and (4.42) and using tlie boundary 
conditions 


P(x = 0, 

D^) 

II 

hJ 

0 

+ ^o' 


Y(x = 0, 

D^) 

II 

< 

0 

+ V ’ 

0 


V(x = L, 

DOJ) 

II 

+ Y’ 

li 


P(x = L, 

D^) 



(4.46) 


the constants and are obtained as 


<^1 = 


T + T ' 
o o 


(4.47) 


H, 


P + P ’ 
o o 

^n04 


G V + Y ’ 
n03 o o 


"n04 




(4.48) 
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Substituting ( 4 . 48 ) in Equations ( 4 . 42 ) and ( 4 . 45 ),- 
the following equations for the velocity and pressure res- 
ponse are obtained. 


Y 


Y + Y' 
o q 


^nx1 


[ 


P +P’ 
o o 


"n04 


C Y +Y' 
n03 0 0 


"n04 E, 


1 


°nx2 ^^*^5) 


Y +Y» 
o o 

E. 


G -z + 
nx5 


[ 


‘^n04 


^n05 

'n04 


Y +Y' 
o o 

E„ 


1 


"'nx4 


( 4 . 50 ) 


Then from Equations (4.49), ( 4 . 50 ) and (4*47), the 
following equations for the velocity and pressure response 
are obtained. 


V^'o 


V h-Y’ C 

o o nL3 

P +P ’ 

00 2 


f — ^ 

L ^n04 


Y +Y' C 
o o n03 


^0-^^; ^n04 ^2 ^ 


(4.51) 


Y +Y' 
o o 


E^ 


P +P » 
o o 

(V +Y‘ ) C 
^ o o^ n04 


•'n03 


"n04 E, 


-I- 


nl 2 
(4.52) 


Dividing Equation ( 4 . 5 I) by Equation (4.52) and rearrang- 
ing the teims yields the following equation 
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P +P> 
o o 

V +V’ 
o o 


Y^+Y'i |G 

nL3 "ni4 ^n03 , 

^n03^m2 


^2 ^2 ^n04 J 

1 ^2=n04 




Px+P*T 


"nP4 

^n04 


(4.55) 


Mere (P^ + P*^)/(Yj^ + Y*^) is the load impedence. 
ii ; method of estimating the load impedence for a load system 
with mercury column oscillating in a vertical U-tube is given 
in Appendix E. The frequency response of the system including 
the effect of longitudinal waves in the tube can be estimated 
from Equations (4.52) and (4.53) with the appropriate system 
parameters. 



CHiPTER Y 


THE EEEEGT OE THE TRMSYERSE WAYES IH THE TUBE 
OH THE DYH1.UIC EESPOIISE OE THE LIQUID 
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The Effect of the Transverse ITaves in the Tube on the 
Dynamic Response of the Liquid: 


given by 


The tube momentum equation in the radial direction is 
18 


, p 

pi (j 

at2 


(A+ 2G) r ir 


i) + (A+ C-) 
r 8r6x 


0^v 

+ G — % + damping 




(5.1) 


The detailed study of Equation (5.1) is given in Appendix 
A for the frequency response anal5^sis. The resulting wave equa- 
tion is linearized by neglecting the "terms of smaller order of 

magnitude by selecting the solution of suitable form and the 

19 

damping effect is included by adding visco-elastic damping . 
The resulting linearized damped wave equation is given by 



(5.2) 


Eor the frequency response analysis, the radial dis- 
turbances are assumed to be of the form 

y(x,t) = y^ e^^ (5.3) 

Then Equation (5.2) is transformed to the ordinary 
differential equation 


9 

d'^y. 


X 


ax" 


+ 0 " 


= ° 


(5.4) 
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whe re 





Assuming the tube is fixed rigidlj^ at both ends, the 
end conditions are given by 

at X = 0 ^x ^ ° 

X = L ° ^5.5) 

Using these boundary conditions the solution of the 
differential equation (5.4) is obtained as 

y^ = EH^sin^ (5.6) 

Uow considering only the excitation caused by the 
fluctuations in the internal pressure, the distribution of the 
above forces on the tube wall can be estimated as given below: 

The internal pressure for viscous flow is given by 
Equation (3.13) as 

Px = E sinh(ir x)+ cosh(ir x^J (5.7) 


Let dx be an elemental length of the tube at any point 

x^. Neglecting the forces acting at all other points except 

at X = x„ 

0 

forcing f inaction = 6(x-x ) Etc R. dx P (5.8) 

Multiplying Equation (5.2) by the area of cross section 
of the tube material and the elemental length dx and, equating 
to the forcing function (5.8) and then rearranging the terms 
yields 
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f ’ 72 
5t^ 


2^- (K + f^) ^ 

OX 


6 ( x-x^ ) 2x Rj_p 


(5.9) 


Taking Laplace transform of Equation (5.9) yieldt 


f' y 




(K+ #1) s 

dx 


2xR. 6(x-x^)p 

^ (5.10) 


By substituting s = jwfor the steady state response 
the above equation reduces to the form 


2 '^^^x , 4 N 2 xR.6(x-x )P 

- G 1 - (i^+ l^j) joJ 1 = i 2L- (5.11) 

^ dx"^ ^ dx^ A 


Combining Equations (5.6) and (5-11) and. rearranging, 
the following equation is' obtained 

o o o 0 o •nTT'V" 2xR-o(x— X )B 

2 H^(-j3»co^+ sin ^ ^ 2 — S (5.12) 


where a^i = G + 0u3( + ^V|) 

Uow by the theory of Eourier series and by the property 
of the dirac delta function is obtained from Equation (5.12) as 


LA + aJ^n^%^/'L^) 


( 5 . 13 ) 


Substituting the expression for in Equation (5.6) the 
radial amplitude of the tube due to the forcing function at any 
point Xq is given by 



(5.14) 


y. 


XX 


where C. 


= P 

3 X, 


sin nux/l sin nixx^/L 
(-/co^+ 


L-i; 


'Transverse oscillations of the tube due to the forcing 
function acting on the entire length of tube is given hy 


y 


X 


L 

^xx 

0 ^o 


dx 

o 


(5.15) 


From Equations (5.14), (5.7) and (5.15) is given by 

.X. "** 


y. 


where C, 


O4 J ¥ 

43 11^ Iq (KRj^) 

P K 


(5.16) 


nrH.(1-C-1)“ cosh (f £))-S,(-1)“ slnlrf;!,)] 

t — 1 (5.17) 

/ ^2 ,y,2 2 /t-2\ /■ , 2, 2 2 2/t-2\ 

(.y +n Ti /L j C-yco + a rn n 71 /I } 


These details are given in Appendix h. 


To predict the radial component of the velocity pertur- 
bations in the liquid, the momentum equation in the radial 
direction, the equation (2.3), is aonsidered. This equation is 
simplified by neglecting the convective acceleration ierms, . 
axial gradients in radial velocity and the radial pressure 
gradients, because these terms are of the second or smaller 
order of magnitude as shovm in Appendix A. The resulting sim- 
plified radial momentum equation is given by the equation. 
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5u’ 

9t “ ^ 


d‘^-a' 4. 1 Su* u* 

8r “ 2 j “ 


S' 8r 


0 


(5.18) 


lov/, if the radial velocity perturbations in the fluid 
are assumed to be of the form 

u' (x, r, t) = u'^ (x,r) e^^ 

Equation (5.19) can be transformed to the form 


4- 


1 


- (4 + ‘): f ) =0 


3r" 


r 8r 






(5.19) 


’®iich is the modified Bessel's equation of order one 
and, for finite solution at the centre of the tube, the solu- 
tion of Equation (5.19) can be written as 

(5.20) 


= g(x, (Kr) 


where K = B 



2^ 

V 


Assuming that the tube remains in contact with the liquid 
the boundary conditions can be written as 


u'^ (x,r=Rj_) = ocoy^ 


n'^ ( r = 0) = 0 


( 5 . 21 ) 

With the above boundary conditions Equation (5.20) reduces 


to the form 


u' 


I^(Kr) 3 ay. 
I^(K%) 


X 


( 5 . 22 ) 


Erom Equation (5.2) it can be seen that the magnitude of the; 
radial expansion is not negligible if the pipe line has relatively 
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low elastic properties. By neglecting the compressihility of 
the liquid^ Equation (2.1) can he written as (Ref .Appendix A) 


6u*^ 

6r 




(5.23) 


Substituting Equation (5.16) in Equation (5.22) and 
then Equation (5.22) in Equation (5.23) the axial flow velocity 
perturbations at any point x is obtained by integrating the 
resulting equation w.r.t.x from 0 to x. 

b 

^’x-'^'x^^’^) = ^5 ^o ^ IT ^ -'•) ^5.24) 

Further it is assumed that t* (0,r) = 0 

where 

D^K G,1 

0 = 1 

^ m I^(KR ) 


From Equation (3.5) the fluid momentum equation for 
the axial velocity perturbations is given by 


9r^ 



9v’ 

X 

Sr 


V 


H dp’ 

1 ^ 

dx 


(5.25) 


By .going through similar steps as used in Equations 
(3.6), (3.7) and (3.8) the solution of Equation (5.25) is 
given by 

1 dp’„ 

v’^ = b(x, I^(Kr) - ^ (5.26) 


and the boundary conditions are 
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3."tl 


r = 0 


6v» 


or 


0 


v’ (x, r = 0) is given by Equation (5.24) 

Substituting above boundary conditions in Equation 
(5.26), the following Equation is obtained 

t (cos^-1) 


dp V 


X 


dx 


3Co>f-[ b(x, 30;) - 0^ S 




n 


(5.27) 


Then Equation (5.26) reduces to the form 


t^(cos _ 1) 


v‘^ = b(x, 3^) [lo(Kr) - 1 J S 


n 


(5.28) 


Average perturbation velocity is obtained by the equation 
_ r p -j \ (cos -1) 

= b(x, 3^) L (KR^) -1 J + E 

(5.29) 


Equation (3.18) can be written as 

dp’ -13 d^’ 

_x _ jB X 


dx 


^ dx^ 


(5.30) 


Substituting Equations (5.27) and (5.29) in Equation 

(5.30) and rearranging the terms yield 


d‘^b(x,.icu) 

dx^ 


2-2 0)2 


Y;2 Ux,M= 05 y-/ 2(14- - 


^^2 

o 


n 

n 


) , cos 


nitx 


^ „ 2 ^5 ^ n 


(5.31) 
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where 


r2 2 
12 ^ 

o 


r 2 " = 




Particular integral of the above nonhomo gene ous ordi- 
nary differential equation is given by 

h 

, = _ n_ r^‘- 5! 

P 


b = - C3 2 


cl? { 2 ^ 


b 


O ^ is cos ^ 
n 1 






(5.32) 


Hence the general solution of the differential equation 
( 5 . 31 ) is given by 


b = G-g cosh(r^x)+ Hg sinh (l^jx)+ b^ (5.33) 

If the average velocity perturbation is assumed to be 
equal to zero at x = 0, then from Equation (5.29) 

b (0, 3^) = 0 (5.34) 

From Equation (5.33) 


b(0, Dw) = % + (O’ 3c4 = 0 

or Gg = -bp ( 0 , jco) (5.35) 

Substituting Equation (5.32) in Equation (5.35) and 
rearranging the terms yield 


G. 


% '2 


4 


n b 


n 

2 ;^ 


n^it 


(5.36) 


( -iT ) 



where 




Kow hy substituting Equations (5.36), (5.32) and (5.33) 
in Equation (5.29) and rearranging the terms yield 


= Or cosh X S —5 5r-o — o — + sinhff^x)- 

^ ^ (y:,^ + n^it^/L^) 7 ^ 1 ^ 


cos mzx/jj 

°22 " n ( r/ . n2.2Ab 
where = ^ ; O^g = 0^ ^ 


— C.T' 'z 2 


33 ^ n 


(5.37) 


G35 = (1 + V ( Cq2)) 


For the pressure perturbations at any point x Equation 
(3.17) can be written as 


p« =im 


(5.38) 


Substituting Equation (5.37) in Equation (5.38), pressure 


perturbations is given by 


X 03 


Gg sinh (f^ S 




P f. 

®7 oj 




^ b^ sin nicx/L 
22 a3L ( r^2^^2^2/^2) 


(5.39) 


Now using the boundary conditions 


at X = L v’^ ” P’x ” 
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in Equations (5.37) and (5.39), the following equation is 
obtained 




V 


_ “nr ^6 (fi^) ^ 


( f^^+n^n^/L^) 
n h 


+ H: 


dBL 


7 oO 


cos 


:h[r:i3: 


/[ Og oosh (f,L) Z ^ r ^J„2/x,2) ^ ‘*7 

c- 1 )“ b. 


C 22 2 


- c^, s 


n 


( n^Ti^L^) 


33 " n 


(5.40) 


Further simplification yields 


H 


7 


r 3B t ,, X B't n ^-n 

l(_^ Cg sinh (f^L)-C^ ^ cosh ({,!)) E 

L (. +n It /L ) 


- G 


P'L 


(-1)^ B 


n 


22 Y« 


+ c,, ^frr^ S ^ 


B ' -rZ . . 2_2/x2n ^ 


n 


( + n^icVL'^) 


r p 1 jB f 

sinh (f^l)- cosh 


Y 


00 




(5.41) 


The frequency response of the velocity together with the 
effect of the transverse waves along the pipe line can be obtained 
by adding Equations (3.22) and (5.37). 


Y = V + V' 

X X 


(5.42) 


By substituting the constant in Equation (5.37) and then 
substituting Equations (5.37) and (3.22) in Equation (5.42) yield 
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or 


where 


V 


Gl Cg cosh (_r^x)S^ + sinh (- 0 ^ 


P’t 
±£ q 

V» ^6 


P’t, 


'33 V' 


S ^) / C^o -°22 


(- 1 )^ cos rnix/L sinh f Ii 


n 


+ S,^ ^2 


10 


^22 V'^ ^9 S 3 ^10 ^ ^ 

(l_(_1 )^ cosh f l) 

"S2^- --'t'-S3So" 


Kg sinh(T 2 ^ 


(5.43) 


V = C ^ G. + G .H. 
nx5 1 nx 6 1 


(5.44) 


°nx5 ° -°6 =°®*" WtISs + Ojil °22 

( - 1 )^ sinh(r l} CO s (nxx/L^ 


n 


+ Srj + Kg cosh(r X ) 


' 11 x 6 


Gg cosh (^f^x^Sg + G^g sinh (K^x)+ Ggg 
cosh(TL)) 


n 


cos “S 3 ^10 "^^2 


P’ 


-G^ + G, 


S. - G. 


G 


n 1 


'9 "5 ‘ "22 Y£: ^6 ~ "33 V*t, 


'10 
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Cg Sg -C 22 Sg + 


n2 


'10 


Pf. 


CO 


- ^6 (fl^)- °6 =oshC^Y) 


^10 (-^1^)“ 


pf^ 


L ' . ' CO 

n^ (-1)^ sinllf l) 




!in]i(f l) 


n 


n 


-7 


d. 


'n 


^ ^ n2(1-(-1)" oosh(fl,)) _ ^ ^ eo,h(fd _ 


■'8 


^10 = ^ 


n 


n 


1_C-1)^ coslxtri.) . 2 . ■ 2 2 2 /^ 2s/ /■2. 2 

— ^ — d ’ ^n n -rcVL )(' +3i tc 

n 


G^ = (-f'o? + nV/L^) ( y^^+n^irV^^) 


at X = 0 


"n05 


= K. ; 


^n06 " ° 


at X = L 


*^nLL5 ~ ~*^6 ^5 cosh (f^ h)+ C^>l sinh ^22 ^6 ^53 ^7 

+ cosh(f l) 

^iiL 6 ~ ^6 ^8 ^n2 sii!Lh(](|Lj + ^22 ^9~^33^10'^ 

Kg sinh(f l) 

Frequency response of the output pressure including 


the effect of the transverse waves along the tube can be 


1 V 3 
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obtained by adding Equations (3.23) and (5.39) 


P = P + u * 

X X 


(5.45) 


By substituting the constant in Equation (5.39) 
and then ' substituting Equations (3.23) and (5.39) in 
Equation (5.45) yield 


P = G. 


sinh 




oBf. cosh/f.x) 


(_Cg + 0^2 yt^ ^ / C, 

■ ^ n(-1)^ sinh(f l) sin nitx/L 


4. n V 

^ ^22 t3L ^ 


s±nh(f 


DBf.C 




1 CO 


sinh(r^x)Sg + 


33 cosh('f^x) 


(CgS8-022 y.^ Sg + S^q) / G. 


IBx n(l-(-l)'^cosh(rL) 

.Q 1 


0 ,d±ui y 
22 COL ^ 


sin + Cg cosh(j x^ 


(5.46) 


or P 


°nx7 Vs Hi 


(5.47) 


where 


^ Og siDh(lC|X)S5 + 0^^ 3^ cosh(L£) 


n(-1)^ sinh(rE)sin rntx/L 

+ 0^2 ^ E — V + °8 - 
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3®^ 33’^. ,, , 

7o~ °6 ■" "w- °n2 =osli(1^x) 


'22 COL 


n(l-(-l)^ CO sh('^ L) sin riTix/L 

^ 3 3 + Og eoshir X 

n 


X = 0 


Pf^ 


*^n07 GO ^n1 ’ ^n08 


pc 

oJ ^n2 ^8 


X = L 


pr^ 


^nL7 ■ a; ^6 % sinh(lf^l)+ cosli (lf^L)+ Gg sini{f l) 

pc .-n^ 

°nI8 = To °6 ^8 (.<1^)+ °n2 Pl °8 !■) 

{JkJ 

The velocity and pressure response including the effect 
of transverse waves along the tube dfe given by Equations (4,47), ; 

(4.48) and (4o4g) in, terms of above constants as 


P, + P ' ' 

Jj Jj 


V + V’ C T-ry 
o q nl7 


p ^ pt p 4 . pt 

0000 


_L_ °n07 

0n08 ■ I’c'*'!’! “n08^2 


(5.48) 




gptf'n 

<*^'oKo8 


°n08 *^2 


Onl.6 (5.49) 


3„+3 ' „ 

o o 

¥ +¥^ 
o 0 




'^nlS ^n07 V , ^n07 ^nI6 ^nl5 
^2 "n08 ^ ^2 '=n08 " ^2 


^nL6 _ L 

^n08 ^1 '^^'l ^n08 


(5.50) 
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liWiere ^ j,) is the load impedence and 

a method of estimating the load impedence is given in Appen- 
dix E. The frequency response of the system including the 
effect of transverse waves can he estimated from Equations 
(5.48) and (5-49) with the appropriate system parameters. 



CHAPTER YI 


CORrPUTEP RESULTS ME LISCUSSIOU 
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6«1 Gom-puter Results : 

Tile physical characteristics given in Table 6.1 are 
used for the computation of theoretical results. 


Tygon tube 


Table 6 . 1 


Load 


Water 


length L ’ . _ 


Inside diameter, 2 

0.1875 inch. 

Outside diameter, 2 R 

0 

0.3125 inch. 

Effective bulk modulus B 

110.1 Ib/in^ 

Density ’ 

1.3 X 62.4 Ib/st 

Poission’s Ratio <T' 

0.45 

Area of U-tube, Aj^ 

0.03267 in^ 

Mass of mercury column, lOj^ 

0.3166 lb 

Spring constant due to gravity, 

per unit area, Kg 

0.98 Ib/in^ 

Density, ■ 

62.4 Ib/ft^ 

Isentropic phase velocity c^ 

90.4 ft/sec 

Kinematic viscosity 

0.00155 in^/sec 


Computer programs for solving the pertinent equations, 
including the programs for the modified Bessel’s function and 
its asymptotic form, were written for the IBM 7044 digital 
computer. 
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To calculate the Telocity and pressure response with 
and without the tube effects, Equations (3.29), (3.30), (4.51), 
(4.52), (5.48) and (5.49) are used. Telocity and pressure 
response without the tube effect were calculated to Icnow the 
magnitude of the effects of the tube on the dynamic flow of 
the liquid. Calculations were made for different Talues of 
the damping constant of the tube material ( and the 

Tiscous damping of the fluid. Knoy/ing the effective bulk modu- 
lus^ modulus of elasticity of the tube material is calculated 
from Equation (2.10). A method of estimating the load impedence 
when the load consists of only an inertia with the restoring 
force obtained from the gravitational field is shown in Appen- 
dix E and the resulting load impedence is given by the equation 





Vp 1 

^ J 


( 6 . 1 ) 


^it is furiher assumed that 



As the experimental results can only be correlated to the 
theoretical results by suitably adjusting some parameters like 
the visco-elastic damping of the tube material, doulomb friction, 
viscous damping etc. the theoretical results are to be calcu- 
lated farr several values of these parameters. For some selected 
parameters the frequency response of the liquid with and without 
the tube effects are given in Figures 6.1 to 6.14. Somtf Figures 
only the effect of the tube on the dynamic response of the 


liquid 
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6*2 G-eneral Discus sion ; 

for the theoretical analysis of the dynamic response of 
liquid flow in pipe lines, linear forms of the basic equation 
of motion have been considered. If the liquid flow/ in smooth 
and uniform circular tube is considered, it is reasonable to 
assume that the fluid flow exhibits axial symmetiy. As dis- 
cussed in Appendix A, the fluid flow equations are said to be 
linear when the disturbances are small in which case the con- 
vective acceleration t erms and the terms consisting of the 
bulk viscosity reduce to extremely small order of magnitude. 

It has also been shoTO in this appendix that for small dia- 
meter tubes the radial flow velocity and the radial pressure 
gradients are of very small order of magnitude. Ihe radial 
momentum equation has been used to include the effect of the 
transverse waves along a soft elastic tube on the dynamic res- 
ponse of the liquid. The energy equation: has not been taken 
into consideration in the entire analysis because of the negli- 
gible changes in the temperature and, hence, in the internal 
energy of liquids due to their high specific heats, low com- 
pressibility and negligible heat transfer to the pipe line system 
when under the action of viscous dissipation. 

For small disturbances in flow conditions, it is reason- 
able to assume' that the equation of state represented by Equa- 
tion ( 2 . 4 ) is isentropic since the dissipation effects in the 
liquid are very small, the specific heats of the liquid and 
the tube material are very high and the energy changes ■' 
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are quite small. Equation (2.4) is a linear form of tlie equa- 
tion for tLe state and the isentropic phase Telocity is given 
hy the linear Equation (2.9). 

It can he seen from Appendix A the assumptions that the 
variations of inside radius and the radial, velocities are veiy 
small are not entirely Justified in the case of fluid flows 
in soft elastic tubes. Though these quantities are resonably 
small they could still introduce some error in the analysis. 

But by indluding the effects in the analysis, the system would 
become quite complicated thus making an analytical study almost 
impossible. Also, in this appendix the fozm of solutions (A. 21a) 
and (A. 21b) make it possible to reduce the equations of motion 
for the tube material to a form suitable for the theoretical 
analysis. Once this is accomplished the terms involving the 
radial gradients vanish and for further analysis it is assumed 
that the form of solution given by Equations (4.3) and (5.3) 
are best suited since the variations in the tube radius are 
quite small . 
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6 *5 Discussion of Results : 

Figures (6,1) to (6.14) show the theoretical c omputed 
results from Equations (3.29), (3.30), (4.51), (4-52), (5.48) 
and (5.49). These figures show that the effect of the longi- 
tudinal waves on the dynamic response of the liquid are very 
small^ compared to the effect of the transverse waves along the 
tube. This is understandable as the shear drag oscillations 
on the tube walls in the case of longitudinal ?/ave are small 
because of the low viscosity of the liquid. 

Prom the above figures it can be observed that the effect of' 
tube. ' reduces as the visco-elastic damping constant increases for 
a particular value of the viscosity of the liquid. It is due 
to the fact that, as the damping constant increases a greater 
portion of the forces on the tube walls is dissipated and hence 
the reduction of the dynamic effects of the tube on the liquid 
flow. It can also be observed that the tube is showing a maximum 
effect around the resonant frequencies of the liquid flow, which 
may be due to the maximum e ffe ct of the forcing function acting 
on the tube at these frequencies. The effect of the tube increase^ 
for a particular value of visco-elastic damping^with the increase 
of viscosity of the liquid, which is due to the increase in the 
effect of the forcing function acting on the tube walls with the 
increase of the viscosity of the liquid. Around the resonant fre- 
quencies of the tube, the effect of the waves along the tube 
does not show considerable effects on the dynamic response of the 
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liquid. This may he, because of the low value of the forcing 
function at the resonant frequencies of the tube owing to the 
dynamics of the main flow. For low values of the visco-elastic 
damping constant, the effect of the transverse waves are not 
shovm. At these values the effect of the tube turn out to be 
very large, hence for these values this analysis is not valid. 
The effect of the tube reduces with the increase of the fre- 


quency 
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STUDY OF THE ODDER OF MIG-IITUDE OF TEEI»IS 


78 


A, 1 A Study of Pluid glow Equations : 

Oscillatory flow of slightly compressible liquid in 
long, elastic, thick walled tubes ?7ith small internal dia- 
meters, is considered for the study of the order of magni- 
tude of terms. Ihe following assumptions can be made for 
this type of system. 



« 1 


V 

o 

c 


<< 1 



<< 1 





<< 1 


(A.1) 


For an elastic soft pipe line system the flow condi- 
tions can be given as follows. 


v^ (the flow velocity) < 3.0 ft/sec. 

I ^length of the tube) f;: 10 ft 
Pw (the inside radius) Ti 0.01 ft 
^ (the resonant frequency of 

oscillations in x direction) 'X. '^*'5 ^ 2% rad/sec. 
c (the phase velocity in liquids) ~ 100 ft/sec. 

Pro ro these values, the order magnitude of the follow- 
ing quantities can be approximately estimated. 
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\ ^ 0.01 
1 10 


.001 


Y 

0 

C 


-X 

100 


.03 


1.3 X 


211 X 0.01 
100 




0.001 


Kie following nondimensional variables are assumed 


for nondimensionalising 

the continui' 

ty and 

momentum 

equations 

Y 

= 


( v’ ) 

^ '^max 

0 

1, 

hvg 

0 

0, 

U 

== 

ev'- 

(u’ ) 

^ '^max 

Q 

1, 

("'hvg 

0 

0, 

P 


Pi 

o 

> 

o 

o 

^max 

0 

1, 


0 

0, 



?oi” ’ 







X 


1 X* , 

^max 

0 

1, 




r 

= 

RirS 

^max 

0 

1, 




t 

= 

1 1- ‘ 

CO ^ ’ 

t’ varies continuously, but the 




duration 

of 

each 

cycle is 

of 

the 


order '[/uJ , that is of period 
1/2 % cc) (A.2) 

vtiere the primes indicate nondimensional quantities 
and the subscript 0 indicates a characteristic constant 
value. In tubes of small internal diameter, the assumptions 
that the radial velocity component is of very small order 
of magnitude as compared to the axial flow velocity is justi- 
fied j particularly so in rigid tubes. The nondim.ensional 
form of the pressure perturbations is assumed on the basis 
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of Joukowsky formula for pressure surges in wave trans- 
mission, where the pressure surge 2\p is given "by 


AP = fcAv 

for a velocity disturbance of ixv. 

Prom Equation (2.1) continuity equation is given by 


1 

f 





Hz + ii |£ + z 

0z ^ _f 6x 


= 0 


(i-.5) 


low substituting the relationships in the set of 
Equations (A. 2) in Equation (A.3) and by multiplying the 
resulting equation throughout by b/v^ yields the following 
equation 

oOL dp * . 16 r 5u* , u’ 1 . 8 v* ^0 v' Bo* 


V T 

0 L 

= R. 


u’ 


Y 6r 





(A.4) 


By considering the terms of the order of unity and 
neglecting the terms of the order of 6 the continuity 
equation in this form with the corresponding terms in 
original form can be v/ritten as 


1_ 8^ ^ 

f "5^ ox 


3u u _ 
8r r 


0 


(1.5) 


Prom Equation (2.2) the momentum equation in the 
axial, direction is given by 


9v , ... Sv . „ iy ^ j! ^ 


0t 8x 8r 




6x 




8^V 


8^v 1 By 

sj 8r2 


0 


( 1 . 6 ) 
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Substituting the relations in (A. 2) and by multiply- 
ing the resulting equation throughout by 1/v^c, yields the 
following equation. 


8v' 
c FF 



V 

Trt 8V’ 


V 

0 , 

av' 

4. 

C 

8x’ 

® 1: 

— u’ 

C 

■ 

a^v' 

+ ^ 


L 1 


ax’ 2 

^ 2 
cRF 

ar*^ 

0 

' 8r’ 


1_ 

f’ 


6p * 

ox’ 


= 0 


(A.7) 


Since the quantities R^/1 and are of the order 

of £ in magnitude and hence, for inertia, forces balanced 
by viscous and pressure forces 


JL a e 

c% “ 


By neglecting the terms of the order of C the momentum 
equation in the axial direction in this form with the corres- 
ponding terms in original form can be written as 


+ 1 Sr 

6t ^ Sx 


- i; 


8^v 1_ 

ar2 r 6r 


1 = 


0 


( 1 . 8 ) 


From Equation (2.3) the momentum equation in the 
radial direction is given by 


j. tr Su , „ 8u 1 ^ 


V 


9^u 
L 9r^ 


+ 


a^u 

ax2 


X ^ 

r Sr 


0 


u 


(1.9) 


Substituting the relations in (A. 2) in Equation (A. 9) 

and , multiplying the resulting equation throughout by 

B/v c , the following equation is obtained 
/ 0 
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£toL dll' 
c 6 1 ’ 


V 6 

+ _2_ y! 

c 


Ml 4. p2 2o ^ 4. Il _L Ml 

ax' . o dr' cr’ 


6L 8 ^u ' €li 1 dll' 

cR^^ 0r'^ cR^^ 


Le u' ^ e e^u* i 

T2 r;2 ^ — 2 I = 0 


(A.10) 


A comparison of Equations (A. 7) and (A.10) shows 
that the terms in the radial momentum equation are of small 
order of magnitude, of order 6 or smaller in magnitude, and 
hence, in the case of a laminar fluid flow in a very rigid 
pipe line, the radial momentum equation can be neglected. 
Prom Equation (A, 10) ^/f ' 6p/8r is of the order S. 

In case of the fluid flow in very soft tubes, v/here 
the radial lelocity and radial momentums may be of consider- 
able magnitude, the simplified radial momentum equation can 
be obtained by considering the terms of the order of mag- 
nitude 6. The radial momentum equation in this form with 
the corresponding terms written in original form is given as 


M 4- 1 ££ 

at ^ f dr 


V 


[ 


d^u 1 u 

arS ^ ^2 


1 = 


0 


(A. 11) 


The radial pressure gradient is of the order of mag- 
nitude 6 and much smaller at still lo?/ frequencies. The 
radial oscillations are governed by a wave velocity which is 


determined by the fluid compressibility and density. These 
facts justify ' the assumption - that the radial pressure 
gradient can be neglected. With this assumption^ the 
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simplified radial momentum equation for liquid flows in an 
elastically soft tube can be written as 



(A.12) 


The simplified form of the continuity equation (A. 5) 

is again considered for further simplification. In the case of 

the liquid flow in a long pipe line v/ith comparatively small 

inside diameter, if the flow is considered as one dimensional, 

the average axial flow can be considered for the analysis of 

the dynamic flow. By taking a cross sectional average of 

Equation (A.5), which can be done by multiplying this equation 

throughout by 2mrdr, integrating betv/een the limits 0 and 

2 

and dividing by the cross sectional 'area m , Equation 
(Ao5) becomes 


2« j- Br eu^u-j j 

' cx or r I 


■jxR^ 


2 0 


71 R,- 


R. (x,t) 

/ 4 Is r d r - 


1 


0 




9t 


11 R- 


e(ixR^‘^) 

at 


(A.13) 


If the values of the average flow velocity and the- 
average density are defined by the equation 


V 


— - / ^ V r d r 

IX R^^ 0 



2ix 

ixRj_^ 


R. 

f ^ f r d r 
0 
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vfriere 


Now 


axLd by 


Equation (A. 13) can be written as 


Q- v(R^,x,t) 8R 

w - “ 


R. 


R.(x,t) 

1 , 2% r 1 8u , u I , 

^ + — / I— +-irdr 


nR, " 0 

1 


JL_ ^ 1 ^ ^ 1 

3t ~ j’ 9t 






6r r 

f(R.,x,t) eR 

6' 


^i 


(A. 14) 


/ ^ 
0 


. 

+ 


- 1 

r i 


n 

rdr = / ""rPdr + u (R^Sx,t) R. 


R- 


r ''i Su , 

- / r ^ dr 


0 


Or 


u(R^,x, t) R£ 


CA.15) 


By substituting Equation (A. 15) into Equation (A, 14) 
assuming 


±_ 

Y 

O 







V 


0 


e 


2 


o 

Equation (Ao14) can 



(A. 16) 

be simplified to the form 

OCitR^^) 
ot 


0 


(A»17) 


ct-lH* 
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Equation (A. 17) can also be written in the form 


6 V 


1 ^ 
B at 


0 


( 1 . 18 ) 


where 


1 

B 




8 >— 



The assumption that the radial velocity and the 
variations in inside radius are small, is partly justified 
from the calculation of these quantitieso 


86 


A, 2 A Study of the Equations for, the Elastic Tubes : 


Tube momentum equation for the longitudinal and 
radial oscillation are given by 


2„ a2_ r-a2. 


r = ( A + 2G) ^ + 


at' 


8x 


(A+ G) [ 


i_2_ + 1 

6x6z 


r 6z 


] 


G 


[ 


9^z . 1 az 


8r 


2 5 ^ 


(1.19) 


6t^ 


(M 2S) + + CA+ 0) 


a rax 


9 

8x 


(A. 20) 


The solutions for the oscillatory motion of the tube 
wall can be assumed as 

_ oia(x - c-]t) (A.21a) 

Z o’ 

y = ^ gia(x - cgt) (A-21b) 

v/here and constants, 

c^ is the longitudinal wave transmission Telocity 
Co is the transverse wave transmission Telocity 

Substituting Equations (A. 21) into Equations (A. 19) and 
(A. 20) and separating the oscillatory components, yields 
the following equations 
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( -a^ = (X + 2G) (-a^) 

, I 


(A.22) 
(A. 23) 


Prom Equations (A.22) and (A. 23) 


. = i'biSS. 

■1 T 


f’ 

Equations (A.22) and (A.23) with the terms ¥/ritten 
in the original :form without damping can he written as 


* ^ = (X + 2G) ^ 


(A. 24) 


at' 


8x" 


8t^ 


(A.25) 


IPPEIIDIX B 


SOilE DE2AILE13 STEPS OE CHAPTER IT 
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® • "I Details of the Steps Between Equations (4.19) an.d (4.20/ 


Equations (4*19) and (4.18) are repeated belo?/ 
1 


X 


f dxo 


0 ^^0 


(B.1) 


XX, 


, , gin ^ sin 

20. h (x^, oco) D I 

‘ ° 2 (B.2) 




Substituting Equation (B.2) in Equation (B.1) yields 

mtx 


20 


„sin 


1 


mix. 


7~;“?T-rr2772; i 3")=i“-Tr *^0 

(-^*ar + an'3i: /I ) 0 


(B.3) 


Erom Equation (3.20) B(x^, juj) is given by 

]i(x^, 3cO) = 0^ cosh (ijx^ + slnh ({ x^^ (B.4) 


nnx 


nixx. 


/ h(x^, 30 j) sin — dx^ = ^ cosii(\ x^ sin— jp- 


dx. 


0 


1 BTCX 

+ / slnh(i^x^ sin — ^ dx^ (B.5) 


The following integrals are taken from the standard 
integral tables'" . 


/ cosh(fxQ^sin — — dx^ = 


nnx 


n7!:3 


fsinh(irx^^sin ^ cosl^x^ cos— ^ 


/2 . ^2^2 /-r2 

1 + n Ti: /L 


(B.6) 


/ siiifc-^jT sin -^—2. 


fcosli(^ 


Xq ) sin 




(f^ + n?x^/L^) 


(3.7) 


Now by substituting Equations (B.6) and (B.7) in 
Equation (B.5) and taking limits gives the follov/ing equation. 


L nxx G (1_(_1)^ cosh(f 

f hCx , joJ*) sin --p-- dx ^ 5— 5- 

0 ° 1 o L ^2 ^ 


(B.8) 

Substituting Equation (B.8) in Equation (B.3) yieldf 


^‘sir.hCf 


n -u ■ 

\ TT 


(B.9) 


where 


20 7c 

= -IT = 


.[ G.,(1-(-l)^cosh{T^H|(-1)^ sinh(fl^] 


(B.10) 
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® ' 2 Detail s of the Step s Betw een Equations - and _( 4-,. JCl 

Equations (4.28), (4.26) and (4.29) are repeated belor; 

(B.11) 


<ll 

M _ 

= Kg g(x 

qco) + 


"d'x 

= D 

[ g (x, 

3^) I 


- 

<|l 


dx 

00 

dx 



(B.13) 


Substituting Equations (B.11) and (B.12) in Equation 
(B.13) yields the following equation 

2 

3«J) I„{K 1 = ^1 

(B.H) 

Now substituting Equation (B.9) in Equation (B.14) 


and rearranging the terms yield 
g(x,jcD) I^(K R.)co2 


o 


OO^ 


^-2 % 


p 2 F 

^o ^ 


DcJCgZb^ sin^ 


= _ QoXIg E b^ sin 


Further rearrangement yields 

_ f g(x,3u) = ^ s (3^1 n! ) b sin ^ 
dx^ ' 2 L 


(B.15) 


,2 ..,2 


f = 1-if ; p2 ^ 


iqCK Ri2 

K. 


where 



iPPMDIX C 


ESTIMATION OF THE EFFECTIVE BULK 
MODULUS OF LIQUIDS IN ELASTIC LIKES 
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Estimation of the Effective Bulk Modulus of Liquids 
in Elastic Lines 


Due to the radial expansion of the tube caused by the 
internal pressure, the voliuEetric capacity of a given length 
of the pipe line is increased. Hence the volume of the fluid 
that has to be introduced into a given length of the tube to 
raise the internal pressure is the sum of the volumetric 
deformation of the fluid due to its compressibility and the 
internal volumetric deformation of the pipe line due to inter- 
nal pressure. 


Considering a pipe line without end restrictions as a 

thick cylinder with inside radius outside radius R^ and 

the Poisson's ratio of the material or, the radial expansion 

of the inside radius R^ for a pressure variation ^^p is given 

*1 S 

by the equation 


AR- 


E 


R.' 

1 




m 


1/E (JAp) 


[^(1 - cr ) + (1 + 0 ") 

(C.1) 


Considering a unit length of the tube, the change in 
fluid volume per unit length caused by the pressure variation 
of Ap is given by 


AV = 


„ [(K. + AS^)2 - E.2] 
% (2R^ + A R^) AR^ 


or 


AV = 


(C.2) 
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For small pressure variations, is neglected 
compared with 2 R^ and hence the variation of the fluid 
volume is given hy 

AV = 2 Ti R^ AR^ (C.3) 


Substituting for AR^^ from Equation (C,1) Equation 
(C,3) reduces to the form 

m. 

AY = 2 n AP (C.4) 

Erom the equation of state (2.4), the volumetric 
strain in the liqiuid due to small pressure variations can 
be written as 


AlSu _ 4£ _ Ap 

Y^ - - f - - B' 


(C.5) 


where B' is the isentropic liquid bulk modulus. 


Hence the total volumetric strain of the liquid and 
pipe line system, which is also equal to the volume of the 
fluid that has to be introduced into a unit internal volume 
of the pipe line system, can be obtained as 


or 


AY 


t 


Y. 




An 


at 


B* ' Y 

2% R,- m 


^ + 

B' 


1 


E 


itR^' 


or 


AY^ 

~Y 


2m. 


t 




(C.6) 


Ihe effective bulk modulus B of the liquid can be 
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obtained as 


or 




AP 


r -1 + 

L B* ER^ 1 


1 - -1 

B " B‘ 


2m, 


ER^ 


( 0 . 7 ) 


This equation can be used to obtain the effective 
bulk modulus when the variations in the inside radius are 
considered to be quite small compared to the inside radius. 



APPENDIX D 


DEIIIDS OE IHE STEPS BETWEEN EQUATIONS 
(5.15) AND (5.16) 
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Details of the Steps Between Equations (5.15) 

and (5.16) 


Equations (5*14) and (5.15) are repeated below 



sin nTtx/l sin nifrx^/D 

(B.1) 

■3 2 

f ,2 „ 2 2 2/^2^ 

(-j3'oo + ttj n Ti /I ) 

1 

N 

o 

dx„ 

. :_o 

(B.2) 


By substituting (D.1) in (D.2) yields 


sin mt 


s/D 


X 


1 

/ P, 


mtiz 


5 (.y.6^ + a 2 0 ^0 


sin 


dx. 


(Do?) 


Erom Equation (3.23)^P^ is given by 

P =^slh sinli(fx^)+ cosh(5rx^] 

^o ^ 

liien 

1 n-rcx^ ' D / r ^ nnx^ 

/ sin dx = Og [g / slnh(fx^ sin — 

0 o n u U Q 

1 nxx 

+ / cosh(fx^sin *^^0 1 (^*4) 

From the integrals given in (B.6) and (B.7) , Equation 


(D.4) reduces to the form 
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L mix 

/ sin 


0 ^0 




C - 1 ) ^slnl:(f ];^H^ ( 1 - ( - 1 ) ^'■' CO sHf. 

(D*53 


2 ^ ^2 2/t2 

+ n % /L 


Substituting Equation (D. 5) in Equation (D. 3) yields 

(E.6) 


^x = ^4 ^ T 


whe re C 


Cg 7. C 3 


4 


pi 


n 


n 


sinhCr:^+ H^(1-C-1)^ coseCl) 
(-|»‘g 4^+ n^')t,VE^)(f^+n'^TC^/E^) 



iPPMmZ E 


LOAD IMPEDMCE MALYSIS FOR A I/IERGUBY 


GOLUIW IN A VERTICAL U-TUBE 



Load. Irfl'psdence Analysis for a Mercury Coluir :.! i;' a 
Vertical U-TuLe 

The dynamic response of the fluid is nueh .inf: u-encsd 

hy the end conditions of the pipe line. Usually the inlet end is 
without any restrictions and it is connected to a supply source. 
There are several possibilities of loading the output or load 
end o.f the pipe line auch as a closed end, an open end, the 
end with a flow r estriction and an end connected to load system 
v/ith characteristics such as inertia, damping, capacitance due 
to coB'pressibility of the liquid and boulomb friction. Load end 
conditions can easily be defined by introducing theierm load 
im,pedence, which is the ratio of the instantaneous pressure 
to the instantaneous velocity of the fluid at the output end of 
the pipe line. The term load admittance is the inverse of the 
load impedence. Por an open exit the load impedance is zero 
and load admittance is infinite, because the gauge pressure at 
the exit is zero. For a closed exit the load impedence is infi- 
nite and load admittance is zero, because the instantaneous ve- 
locity component at the output end is zero. Srinivas has given 
a method of estimating the load impedence for a load system ?/ith 
mercury column oscillating in a vertical U-tube and is given as 
follows . 

Neglecting the damping effect of mercury in a glass 
tube, the force acting on the mercury column and the load 
pressure, inertia of the column of mercury and the restor- 
ing force due to the gravitational field. Hence the force 
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balance eq,nation can be vnritten as 


For frequency response analysis the time derivative 
and the time integration in Equation (E.1) is replaced by 
2 cO and l/;ia3 . By the above substitution Equation (E,1) 
reduces to the form 

h = ^ [ j CO H- ^ ] Wj, (E.2) 

The ratio of the velocity of the liquid at the pipe 
line and at the load is given by 




A 

_e 


and hence the load impedence of this system is obtained as 




1 

_S p 

e^Aj^ 


] 


CE.3) 



